Casimir manipulations: The orientation dependence of fluctuation-induced forces 
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The Casimir interaction between two objects, or between an object and a plane, depends on their 
relative orientations. We make these angular dependences explicit by considering prolate or oblate 
spheroids. The variation with orientation is calculated exactly at asymptotically large distances for 
the electromagnetic field, and at arbitrary separations for a scalar field. For a spheroid in front 
of a mirror, the leading term is orientation independent, and we find the optimal orientation from 
computations at higher order. 
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Casimir forces, predicted in 1948 [T] and probed in high 
precision experiments over the past decade [H [31 IH [S], 
are particularly important at micro-meter to nano-meter 
length scales. In constructing and manipulating devices 
at these length scales it is important to have an accurate 
understanding of the material, shape and orientation de- 
pendence of these forces. Since these dependencies of- 
ten originate in the inherent many-body character of the 
force, they cannot be reliably obtained by commonly used 
approximations. In this Letter we investigate shape and 
orientation dependence of the Casimir force using a re- 
cently developed method [S] (see also [7]) that makes it 
possible to compute the Casimir interaction for arbitrary 
compact objects based on their scattering matrices. 

A closely related work is the classic paper of Balian & 
Duplantier [5], which provides expressions for the large 
distance asymptotic force between perfectly conducting 
objects of arbitrary shape and orientation. Using our 
method, we generalize these results in a form suitable for 
extension to arbitrary dielectrics and distances. As tan- 
gible examples we then focus on ellipsoids, computing the 
orientation dependent force between two spheroids, and 
between a spheroid and a plane. The latter geometry is 
particularly interesting, as there in no orientation depen- 
dence at leading order for any object (metal or dielectric) 
in front of a mirror. We obtain the preferred orienta- 
tion of a spheroid from a computation at higher order, in 
which T-matrices are evaluated in an expansion around 
a spherical shape. To numerically study a wider range of 
eccentricities and separations, we consider a scalar theory 
(in place of electromagnetism) , for which efficient com- 
putations can be performed using a basis of spheroidal 
harmonics. For the case of Neumann boundary condi- 
tions, we find a transition in the preferred orientation as 
the spheroid approaches the plane. 

We start from a general expression for the Casimir en- 



ergy as an integral over imaginary wave number f6|, of 
Tr ln(l — N), which after expansion in powers of N reads 
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Here N = T^V^'^T'^H'^'^ , where and relate incom- 
ing and scattered electromagnetic (EM) fields, while the 
"translation matrix" U^^ relates the incoming wave at 
one object to the outgoing wave at the other. Natu- 
rally, evaluating energies (and forces) from these formal 
expressions requires judicious choice of basis. 

A natural basis for the field is (vector) spherical har- 
monic multipoles centered at each object, and labeled 
by (/,TO, A); the additional index A G {M, i5} distin- 
guishes between transverse electric (TE) or magnetic 
(TM) modes. Each factor of U decays as e"'*'^, so we can 
obtain an asymptotic series, valid at large separations 
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\ ^ a, the expansion is dominated by the lowest order 
multipoles. Because electromagnetism does not admit 
monopole fluctuations, the leading asymptotic behavior 
comes from p- waves, I = I' = 1. By examining the re- 
sponse of an object to a uniform electric (magnetic) field, 
it is straightforward to relate elements of the 3x3 ma- 
trix "^fmim' C^imlm') elements of the standard elec- 
tric (magnetic) polarizability matrix a {j3). In terms of 
the Cartesian components of the latter (with the z axis 
pointing from one object to the other), the first asymp- 
totic contribution to the energy is 



^7 I KP^xx^xx 



yy yy 



2al,aD (2) 



xy^xy 
yz^yz) + P) 

7 {alj^yy + alyPl, - 2aly/3ly) + (1 ^ 2) 1 , 



30 (al.„a 



2 

xz^xz 



2 



which generalizes the result of Balian and Duplantier [5] 
for perfectly conducting objects at nonzero temperatures. 
For the case of an ellipsoidal object with static electric 
permittivity e and magnetic permeability fj,, the polar- 
izability tensors are diagonal in a basis oriented to its 
principal axes, with elements (for i £ {1,2,3}) 
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where V — 47rrir2r3/3 is the ellipsoid's volume and the 
so-called depolarizing factors are given by 
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in terms of the semi-axis dimensions r^. We can gain 
further insights into these results by focusing on the case 
of spheroids, for which ri = r2 — R and ^3 — L/2. 
Then the depolarizing factors can be expressed in terms 
of elementary functions. 
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where the eccentricity e = y 1 — is real for a prolate 
spheroid {L > 2R) and imaginary for an oblate spheroid 
(L < 2R) . The polarizability tensors for an arbitrary ori- 
entation are then obtained as a = TZ~^a^TZ, where TZ is 
the matrix that rotates the principal axis of the spheroid 
to the Cartesian basis, i.e. 7?.(1,2,3) — > {x,y,z). The 
result provides the leading Casimir energy for spheroids 
of arbitrary size, orientation, and material. Note that 
for rarefied media with e ~ 1, /i ~ 1 the polarizabilities 
are isotropic and proportional to the volume. Hence, to 
leading order in e — 1 the interaction is orientation inde- 
pendent at asymptotically large separations, as we would 
expect, since pairwise summation is valid for e — 1 <C 1. 

As a special case, we present the explicit formula for 
identical perfectly conducting prolate spheroids with the 
length L ^ R. To implement the perfect conductor limit, 
we take e ^ 00. In this limit, the full T-matrix is inde- 
pendent of /i. Because we are doing an expansion around 
zero frequency, however, we must set /i = to suppress 
the magnetic field that is still allowed at zero frequency. 
The orientation of each "needle" relative to the line join- 
ing them (the initial z-axis) is parameterized by the two 
angles {9,^), as depicted in Fig. [l] Then 
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FIG. 1: Orientation of a prolate (cigar-shaped) spheroid: The 
symmetry axis (initially the z-axis) is rotated by 9 about the 
X-axis and then by tf) about the z-axis. For two such spheroids, 
the energy at large distances is give by Eq. The latter is 
depicted at fixed distance d, and for -01 = i/'2, by a contour 
plot as function of the angles 61 , 62 for the x-axis rotations . 
Minima (maxima) are marked by filled (open) dots. 




FIG. 2: As in Fig. [T] for oblate (pancake-shaped) spheroids, 
with a contour plot of energy at large separations. 



We note the following features of this result: 

• The energy is minimized for 9i — 92 — 0, i.e., for two 
needles aligned parallel to their separation vector, as we 
expect since this case is dominated by the fluctuations of 
electric dipoles along the symmetry axis. 

• At almost all orientations the energy scales as L^, and 
vanishes logarithmically slowly as i? — > 0. 

• The dependence vanishes when one needle is orthogonal 
to z (i.e. 01 = 7r/2), while the other is either parallel to 
z {02 — 0) or has an arbitrary 02 but differs by an angle 
7r/2 in its rotation about the z-axis (i.e. ipi — 4^2 — i'/2). 
In these cases the energy comes from the next order term 
in Eq. and takes the form 
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• From Eq. we see that the least favorable configu- 
ration (6*2 — 7r/2) corresponds to two needles orthogonal 
to each other and to the line joining them. 

For perfectly conducting oblate spheroids with R ^ 
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L/2, the orientation of each "pancake" is again described 
by a pair of angles {9, as depicted in Fig.|2] A contour 
plot of the leading angular dependence is also presented 
in the figure; its explicit formula will not be reproduced 
here, but we note the following features: 

• The leading dependence is proportional to R^, and does 
not disappear for any choice of orientations. Further- 
more, this dependence remains even as the thickness of 
the pancake is taken to zero {L — + 0). This is very dif- 
ferent from the case of the needles, where the interaction 
energy vanishes with thickness as ln~^ (L/R). We at- 
tribute this effect to dimensionality: The limiting linear 
needle is only marginally visible to the EM field, while 
the limiting two dimensional plane remains as an opaque 
obstacle. The lack of L dependence is due to the assumed 
perfectly metallic screening. If the dielectric function re- 
mains finite in the static limit, a conventional scaling 
with volumes of the objects is expected. 

• The configuration of minimal energy corresponds to 
two pancakes lying on the same plane (^i = O2 — 7r/2, 
ip ^ 0) and has energy -he {173 / ISn^) R^ /d\ This is 
due to the electric dipole fluctuations in the planes. 

• When the two pancakes are stacked on top of each 
other, the energy is increased to —he{62/9n^)R^/(f. 

• The least favorable configuration is when the pancakes 
lie in perpendicular planes, i.e., 9i — n/2, O2 — 0, with 
an energy -he{n/-iTT^)R^ /(f . 

For an object interacting with a perfectly reflecting 
mirror, we can use its image to construct N — — T^U^'^^, 
where R stands for the reflected object [10 . Alterna- 
tively, we can express the scattering matrix for the mirror 
in a basis of plane waves, and use a translation matrix to 
convert between this basis and the basis for scattering of 
the compact object. At leading order we have 
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which is clearly independent of orientation. A similar ex- 
pression for a mirror and an atom is present in the classic 
work of Casimir and Polder [T], and implicit in Ref. 
for perfect conductors. (The orientation independence of 
energy is special to zero temperature; at nonzero temper- 
atures the general expression in Ref. |8j does contain an- 
gular dependence.) Orientation dependence in this sys- 
tem thus comes from higher multipoles. The next order 
also vanishes, so the leading term is the I = 3 contribu- 
tion to the term linear in N in Eq. ([T|), in which case the 
scattering matrix is not known analytically. 

Instead, we obtained the preferred orientation by con- 
sidering a distorted sphere in which the radius R is de- 
formed to R + Sf{-d, (fi). The function / can be expanded 
into spherical harmonics Yim{'d,^), and spheroidal sym- 
metry can be mimicked by choosing / — 1^20(^1 V'). 
By performing a perturbative expansion of spherical T- 
matrices in S, the leading orientation dependent part of 



the energy is obtained as 
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A prolate spheroid (5 > 0) thus minimizes its energy 
by pointing towards the mirror, while an oblate spheroid 
((5 < 0) prefers to lie in a plane perpendicular to the 
mirror. (We assume that the perturbative results are 
not changed for large distortions.) These configurations 
are also preferred at small distances d, since (at fixed 
distance to the center) the object reorients to minimize 
the closest separation. 
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FIG. 3: The orientation dependence of the Casimir energy 
for a prolate Neumann spheroid with LjR — 4.8 opposite a 
Neumann plane, as a function of the angle 9 between the nor- 
mal to the plane and the spheroid axis, for various values of 
d, the distance from the spheroid center to the plane. The 
asymptotic result of Eq. ( 12 1 is depicted by the top line. For 



d < 8.5R, the spheroid prefers to point toward the plane, 
while at larger distances it spontaneously breaks the rota- 
tional symmetry and aligns parallel to the plane. The larger 
the value of ^0 — 1/e, the larger the distance at which this 
transition occurs. 

Interestingly, the above conclusion is not generally 
true, and as discussed below, we find that there can be 
a transition in preferred orientation as a function of d 
in the simpler case of a scalar field. Then we can ob- 
tain the T-matrices and translation matrices as straight- 
forward generalizations of the spherical case, in which 
spherical harmonics and Bessel functions are replaced 
by spheroidal harmonics and spheroidal radial functions. 
Prolate spheroidal coordinates ^Tj are defined by taking 

as the usual azimuthal angle and 



= a^ie-m-v'), z = aTj^, (10) 



where the interlocal separation is 2a. The surface of the 
spheroid is defined by a constant spheroidal radius ^ = 
^0 = 1/e, with L = 2^oa and R — a^/^o ~ 1- Oblate 
coordinates are obtained by analytic continuation of the 



4 



radial coordinate with ^ ^ — and k ^ ik. We use the 
conventions of Ref . [12j in which the spheroidal functions 
are normalized identically to their spherical analogs. 

The T-matrix is diagonal, with elements which in 
terms of radial spheroidal functions of the first and third 



kinds are Tf^{iK) 



"(inai^o) 



^(3), 



, where the prime in- 



dicates derivative with respect to the second argument. 
The translation matrix converts between the spheroidal 
basis, indexed by I and m, and the plane wave basis, in- 
dexed by the momentum parallel to the plane q±. The 
calculation is performed on the imaginary axis; with q± 
Teal, and gy = iy^n"^ + q^, which is pure imaginary. The 
translation matrix element Uj^^ (ik) at separation L is 

UL''^' = t'i-ir ^y^'""{^Ka, 7^, (j)), where is a 

1 -111 • 1 n I] cos 0—q± cos (ha sin 6 -tt 

spheroidal harmonic and r/ = — — r« ^ Here, 

9 is the angle between the spheroid axis and the normal 
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to the plane and tan < 
is the polar angle of q±^ around the perpendicular to the 
plane. To implement the matrix multiplication needed 
to compute N, we integrate over q^. Normalizing the 
integration measure introduces a factor of 2q\\/{iK). 

The numerical results can be compared to the large 
distance expansion for an object with Neumann bound- 
ary conditions that is placed opposite to a plane (the 
latter with Dirichlet or Neumann boundary conditions). 
The general form of the energy in this case is 
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for a Dirichlet/Neumann (+/— ) plane, where V is the 
volume, and /? is like the static magnetic polarizability 
tensor in Eq. ([s]) for /i = 0. For a spheroid this results in 

fie R'^L 9-An3-nl + (Sng - 1) cos 29 
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where is given by Eq. ([5]). In contrast to the EM en- 
ergy of Eq. (|8| , the leading term of Eq. ( 12 I is orientation 
dependent. For a prolate spheroid with R ^ L we have 
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and for an oblate spheroid with R^ L, we have 
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To leading order in R/L and at fixed separation d, the 
energy for a Neumann plate is minimized at 9 — tt/2 for 
a prolate spheroid and at 9 = for an oblate spheroid. 
This means that both a needle and a pancake prefer to 
be parallel to the plate; in the former case because of 
Pxx and Pyy and in the latter case because of Pzz- For 
a Dirichlet plate, the energy is minimal if the needle or 



the pancake are perpendicular to the plate as would have 
been predicted by pairwise summation. 

In Fig. |3] we present a sample of the full numerical cal- 
culation, using an optimized version pi3j of the spheroidal 
harmonic package of Falloon [T3]. We truncate the infi- 
nite sum at I ~ 2, though we have verified for individual 
cases that summing through I — 4 does not significantly 
change these results. At close separations, the prolate 
spheroid prefers to point toward the plane, as we expect 
from the proximity force approximation. At large sepa- 
rations, we have shown that the configuration parallel to 
the plane has the lowest energy. We thus observe a spon- 
taneous breaking of the symmetry around the perpendic- 
ular to the plane; the separation at which this transition 
occurs varies with the spheroid's eccentricity. 

These examples provide a glimpse into the myriad ef- 
fects of shape and orientation on the Casimir effect. One 
can imagine that with better precision of experiment and 
theory such forces can be employed to manipulate small 
objects at the micro- and nano-meter scales. 
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